Simulation of the thermally induced austenitic phase transition in NiTi nanoparticles 
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The reverse martensitic ( "austenitic" ) transformation upon heating of equiatomic nickel-titanium 
nanoparticles with diameters between 4 and 17 nm is analyzed by means of molecular-dynamics 
simulations with a semi-empirical model potential. After constructing an appropriate order pa- 
rameter to distinguish locally between the monoclinic B19' at low and the cubic B2 structure at 
high temperatures, the process of the phase transition is visualized. This shows a heterogeneous 
nucleation of austenite at the surface of the particles, which propagates to the interior by plane 
sliding, explaining a difference in austenite start and end temperatures. Their absolute values and 
dependence on particle diameter are obtained and related to calculations of the surface induced size 
dependence of the difference in free energy between austenite and martensite. 

PACS numbers: 02.70.Ns, 64.70.Nd, 68.35.Rh 

I. INTRODUCTION 

Shape-memory alloys (SMAs) exhibit the unique effect of "remembering" their original shape after deformation 
and heating ["shape memory effect" (SME)], and moreover, they show "super-" or "pseudoelastic" behavior (SE), 
i.e. recoverable strains of up to S%X These features make them to widely used and studied materials in industry 
and research. The underlying physical effect is the martensitic phase transition (MT), where a crystal lattice trans- 
forms diffusionless between a high-symmetric structure ("austenite") and a low-symmetric structure ("martensite"), 
depending on external temperature or stress. Within the scope of nanoscience and -technology, extensive effort has 
been made in the last decade to explore MT, SME, and SE at the nanoscale (see Waitz et al^ for an actual overview). 
Especially in the case of near-equiatomic nickel-titanium, the most frequently used SMA, many experimental studies 
examine the behavior of different nanostruetures such as particles^^— , thin films^iiS or pillars^i. In most of these 
works, the dependence of system size on the structural transformation temperatures or stresses on the one hand, and 
on resulting morphologies on the other hand is discussed, as well as the question, whether there is a critical diameter, 
below which no MT can occur. 

In order to analyze this kind of questions from a theoretical point of view, molecular dynamics (MD) simulations 
have been performed in the past for different nanostructured alloys and pure metals, which can also show MT. Exam- 
ples are simulations of Cu filmsi^, NiAl wiresi^, and CdSe partielesi^, where stress induced transitions take place, or 
of FeNi particlesi^ and Fe wiresi^ in the case of thermally induced MT. Visualization of the atomistic transformation 
behavior of these systems has shown, that mainly free surfaces act as nucleation sites for the emerging structures, 
suggesting an explanation of decreasing transition temperatures with system size by an increasing surface-to- volume 
fraction. 

This work focuses on the thermally induced reverse martensitic ( "austenitic" ) phase transition in NiTi nanoparti- 
cles, occurring between the monoclinic B19' at low temperatures and the cubic B2 structure at high temperatures. 
MD is applied with a Finnis-Sinelair potentialii, with NiTi parameters determined by Lai and Lii*i^. In a recent 
worki^, the authors of the present paper applied these methods successfully to bulk NiTi, where a strong dependence 
of transition temperatures on Ni concentration has been found and could be explained by energetic calculations. The 
present work extends these studies to equiatomic NiTi nanoparticles, by visualizing the transformation process, and 
by determining for the first time absolute values of austenitic transition temperatures of NiTi as a function of system 
size. As in the case of concentration variation, this dependence can be explained by the change of the energy of the 
system, in this case due to free surfaces. The paper is organized as follows: In SecHH a detailed structural analysis 
is presented and an order parameter is defined, which allows to distinguish between the occurring structures. In 
addition, details of the simulation method are given. The phase transition upon heating of the particles is analyzed 
and discussed in See. lIIII Finally, a conclusion is given in Sec.lIVl 

II. STRUCTURAL ANALYSIS 

As mentioned above, it has been shown recentlyi^, that a slight modification of the cut-off behavior of effective 
hopping integrals and pairwise interactions leads to a stable B19' structure with improved lattice parameters a — 
4.45 A, 6 = 4.03 A, c = 3.00 A, and a = 97.7° of the monoclinic cell at T = K in bulk systems, if compared to 
experimentalSS and ab initio^ values (see Fig.[TJa) for a schematic structure representation). The body-centered 
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orthorhombic structure B33, proposed theoretically by Huang et ali^S to be the lowest energy state of NiTi with a 
shear angle of 107° is not reproduced by the choice of potential parameters used in this work. Until now, B33 could 
not be stabilized experimentally in pure NiTi, neither in bulk nor in nanoscaled systems, which is assumed to be due 
to internal stresses, leading to the monoclinic B19' as the martensitic structure^. Besides B19' at low temperatures, 
the potential used in this work predicts a phase transformation to a cubic structure close to B2 (c = 3.01 A) if the 
system is heated above the austenitic transition temperature Ta = 318 K in aMD simulation under periodic boundary 
conditional^. Since the reverse transition upon cooling produces a related, but different crystal structure, which 
is energetically slightly more favorable than B19', the present work focuses only on the austenitic phase transition 
process. The nanoparticlcs are realized as spheres under free boundary conditions. MD is performed by a velocity- 
Verlet algorithni^^ (timestep Ai = 10~^^ s) with temperature control by a Nose-Hoover thermostat^. 

Since it is not possible to identify the occurrence of the structural phase transition by changes of the shape of 
a simulation box when applying free boundary conditions, an order parameter is necessary to distinguish locally 
between B19' and B2. Therefore, the nearest neighbor environments of these two structures are considered, as they 
are produced by the potential in a bulk simulation. In B19' [Fig.[IJb), left], the dashed lines symbolize the neighbor 
distances, which are elongated about 13.5% relative to the continuous ones due to the monoclinic shear in [100] 
direction. In B2 [Fig.[IJb), right], the formerly two sheared lengths have nearly the same values as the unsheared 

ones [there is a difference in neighbor lengths of ±3% relative to a perfect B2 NiTi (2.62 A), which is denoted by 
the dashed and continuous lines in this case] . By denoting the mean value of the two lengths in [100] direction as di 
and the mean value of the six lengths in the other directions as do, one can construct a quantity x = X {do, di), which 
shaU be equal to 1, if di = df^^' = 2.85 A and do = d^^^' = 2.51 A, and equal to -1, if di = df ^ = do = d^^ ^ 2.62 A. A 
linear ansatz x {do, di) ^ Aoda + Aidi leads to: 

do (d^2 + df"') (dS2 + d,f«') 

This order parameter is capable of distinguishing between B19' and B2 on an atomic scale and of resolving intermediate 
states, but it is restricted to single variant domains with known shear direction, as studied in this work. By looking 
at systems with differently orientated monoclinic variants, one has to search for the two opposite elongated next 
neighbor distances in B19' for each atom first, before applying Eq.[T] 

In order to simulate the austenitic phase transition in equiatomic NiTi nanoparticles, they are set up as spheres 
with diameters d ranging from 4 to 17 nanometers in a perfect B19' single variant structure, corresponding to a 
number of atoms between about 2000 and 196000. In this size range, the single B19' structure in the nanoparticles is 
stable, but since the total energy of the surface atoms in B2 structure is lower than in B19' orientation, this is only 
valid below a critical surface-to-volume ratio corresponding to a minimal size. Below this size, the particle contains 
of B2 even at low temperatures. After an equilibration time of 40000 timesteps, where the imposed temperature is 
reduced stepwise to 1 K, the systems are heated about 2 K every 2000 simulation steps up to 350 K. Heating rates 
in the range of 1 K/ps are usual in MD simulations of solid-solid phase transition a^^'^^ , and can even be realized in 
experiments by laser irradiation and shock-wave loading2L2^. By performing MD, Sandoval and Urbassek showed, 
that there is only a weak dependence of transition temperatures of Fe nano wires on the heating/cooling rate^^. 

III. RESULTS AND DISCUSSION 

Fig. [5] shows the mean x- value of all atoms during the heating process at 4 different system sizes. The smaller 
the particle, the more atoms lie at the surface, leading to the deviation of x in the martensitic state from the bulk 
value x=li since the surface atoms do not "see" a perfect B19' structure around them. At first, x increases and 
can obtain even values above 1. This is due to thermal expansion of the next neighbor distances leading to the 
atomic values (see Eq.[T|), out of which x is calculated. Nevertheless, the structures remain stable until a certain 
temperature is reached, where a transformation to the B2 phase (x = -l) starts. The difference in austenitic start and 
end temperatures (called As and Af) indicates a heterogeneous transformation, starting at certain nucleation sites and 
propagating through the system. Since this difference depends on particle size, a surface triggered transformation can 
be assumed. Therefore it is necessary to use a more local order parameter instead of the global x, where the core atoms 
are not taken into account, to identify the onset of the structural change more precisely. To this end, the nanoparticles 
are subdivided into spherical shells with a thickness of 1 nm, and the mean x in each of these shells is calculated. For 
the particle with d = 8 nm, the result is shown in the inset of Fig. [51 where the different temperatures, at which the 
austenite reaches the inner lying shells, are visible. With this method. As can be determined quantitatively for all 
system sizes, as the temperature, at which a significant drop of the structural order parameter of the outermost shell 
occurs. Af is defined as the temperature, where the x-value of the innermost shell reaches -1. Values x<-l in the 
vicinity oi Af are due to non-equilibrium states slightly deviating from B2. 
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FIG. 1: To distinguish locally between B19' (left half) and B2 (right half): (a) schematic representation of the lattice cells; 
(b) nearest neighbor environments with atoms in shear direction denoted by s and neighbor distances 2.85 A (B19', dashed), 
2.51 A (B19', continuous), 2.69 A (B2, dashed), 2.54 A (B2, continuous); (c) resulting value of the order parameter x ^nd a 
translation to a color code 




FIG. 2: Mean value of the structural order parameter x of ^-U atoms during a heating process of nanoparticles with different 
diameters d. Inset: Mean values of x in spherical shells with a thickness of 1 nm. 
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FIG. 3: Visualization of a transforming nanoparticle with diameter d=10 nm. The xz-plane through the origin is shown at 
temperatures (A) 239 K, (B) 259 K, (C) 275 K, and (D) 289 K. 



TABLE I: Differences in surface energy per atom between B19' and B2 [Ae'''*^'' = e^''*''' (B19') - e^'"'"'' (B2)] for some lattice 
planes (hkl) (in eV). 



Lattice planes (hkl) 




(010), (111), (111) 


-0.06 


(ioo),(iii),(iii) 


-0.12 


(001),(011),(0ll) , (110),(101),(I10) , (101) 


-0.17 - -0.21 



In order to understand the atomistic details of the transformation behavior, a visualization of the system is done 
by relating the atomic % values to colors [Fig.[TJc)]. Fig. [3] shows a central plane of a nanoparticle with diameter 
d=10 nm at four different temperatures between Ag and Af. The austenite transformation starts at the surface 
and propagates through the particle by sliding of certain lattice planes. If the particle is visualized not only in the 
xz-plane, but also in the xy- and yz-planes in the same manner, a martensite polyhedron can be identified, which is 
surrounded by (100), (111), (111), (HI), and (111) lattice planes. This behavior can be understood as follows. The 
atoms at the surface of a spherical nanoparticle are arranged in many differently orientated planes. Each of these 
free surface planes has a certain energy difference per atom between martensite and austenite (Ae). If the particle 
is heated, the atoms at planes with the smallest Ae values will change their structure first, at T = As. TableU shows 
these differences for certain lattice planes, calculated separately by setting up systems in B19' and B2 structure with 
free surfaces of the desired orientation at T = K, equilibrating them in a simulation, and averaging the energies of 
the surface atoms. 

In general, energetic barriers have to be overcome in martensitic transforming materials to enable the first order phase 
transitions. These barriers result from nonchemical energy contributions as elastic strain energies or unlike interfaces, 
which necessarily have to emerge, if a certain crystal structure forms out of another onei. This barrier, which prevents 
a transformation below a certain temperature in the bulk like interior of the particle, is lowered for the atoms lying 
at the interfaces between the already transformed surfaces and the bulk. This leads to a structural change of these 
atoms at a lower temperature compared to the bulk value of A^, what is proceeded plane by plane, until a temperature 
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FIG. 4: Dependence of austenitic transformation temperatures on particle diameter and the total number of atoms. Inset: 
difference between austenite and martensite energy per atom as a function of 1/d. Lines represent linear data fits. 



slightly below Af is reached. Here, the remaining martensite transforms homogeneously to B2, which can also be 
seen in the x-curves of Fig. [51 where the slopes drop nearly vertically before the transformation is completed. Going 
to larger systems, the influence of the surface shrinks, which implies a more homogeneous transformation and a less 
pronounced lattice plane sliding, resulting in a decreasing difference between Ag and Af. 

Plotting the simulation results of As and Af as well as the mean austenitic transition temperature Ta~0-5 
{As + Af) against l/d (see Fig.|4]) leads to the behavior 

Ta - T^""^ - ^ (2) 

by linear regression, with parameters /3 = 536.3 nm-K, and T^""^ = 309.11 K. Thermodynamic calculations of different 
nano particles, which undergo structural phase transitions, have shown, that in these systems, the difference in free 
energy per atom between the two involved structures (A/""™) shows a proportionality with 1/d, too ^^i"^" . Since 
this result is obtained by separating A/""™ into a bulk and a surface part, it can be assumed, that it is the surface 
induced size dependence of A/""™, which determines the 1/d behavior of the transition temperatures. For the NiTi 
nanoparticles, a calculation of A/""™ at T = K, where A/""™ equals the difference in internal energy per atom, has 
been carried out, which confirms the linear dependence on 1/d (see inset in Fig. 2]). The results are obtained by setting 
up different sized spheres in B19' and B2 structure, and calculating the mean energy per atom after equilibration at 
T = K. 

The T^""^ value of Eq.[2]fits well with Ta = 318 K, obtained by periodic boundary conditionsi^. Because in that 
case, the transformation is completely homogeneous, no difference between As and Af could be identified. But since 
real NiTi systems always contain free surfaces in the form of grain boundaries, as well as other lattice defects or 
dislocations, which can act as nucleation sites, one should expect a non-zero A^ ^Af — Ag, what in fact is observed 
experimentally in bulk NiTi with Ni concentration between 50.3% and 51%, lying in the range of 29-45 K^. This 
is in good agreement with the values obtained for nanoparticles bigger than d^vll nm (50000 atoms), for which A^ 
seems to stay within about 37 ± 5 K. Since the absolute experimental value of Ta (331 K for Ni5o.3Ti49.7)^ fits well 
with T^""', it is assumed, that overheating effects due to the strong heating rate of 1 K/ps are negligible. 

Eq.[5] was observed in other nanosized systems and materials showing structural phase transitions, too, such as 
Co particles^, FeNi particlesi^ or Fe nanowiresi^. In general, the proportionality of transformation temperature to 
the inverse diameter in small systems, which undergo first order phase transitions is known as the "thermodynamic 
size effect". It was derived at first by Pawlow in 1909 for the solid-liquid transition in small particles by a pure 
thermodynamic calculation of the gas pressure at a given particle size and temperature^^. 
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IV. CONCLUSION 



In conclusion, an analysis of the austenitic phase transition in spherical NiTi nanoparticles is presented. By 
calculating next neighbor lengths, an order parameter is constructed, which can resolve changes between the low 
temperature structure B19' and the high temperature structure B2 at an atomic scale. It is shown, that there is a 
heterogeneous structural transformation from the outer to the inner lying atoms of the particles upon heating above the 
austenitic start temperature Ag. Visualization of the transformation process shows a sliding of martensite/austenite 
interfaces of certain orientation, triggered by unequal energy differences between B19' and B2 configurations at 
different free surface planes. Transition temperatures as a function of system size are obtained, leading to a T^i oc 
law, which is well established in nano science. This behavior is related to the size dependent energetic difference 
between austenite and martensite, triggered by the surface-to- volume ratio of the nanoparticles. 
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